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, Abstract. I review the status of next-to-leading-order calculations for hadronic 

final states in deeply-inelastic lepton-nucleon scattering. In more detail, I focus on 
calculations of (2-|-l)-jet-type cross sections, describe recent progress in extending 
the perturbative description into the transition region between deeply inelastic 
scattering and photoproduction, and study the validity of the fragmentation 
function picture for one-particle-inclusive cross sections at small and small Xp. 
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' 1. Introduction 

cn . 

, Recent experimental results from HERA at DESY show that the hadronic final state in 

deeply inelastic scattering can be studied with high precision. The results include the 
measurement of the strong coupling constant as{Q^) by means of the (2-|-l)-jet rate 
i?2-i-i = cr2+i/crtot [0, H, H and event shapes a direct determination of the gluon 
density [pi, and the measurement of momentum fraction distributions for charged 
pi ' particles Q . The latter indicates that it may be possible to study scaling violations 

D , of fragmentation functions, the virtuality of the photon being the relevant scale 

' for the fragmentation process. 



In this review I will concentrate on three selected topics: 



j_j . • NLO calculations for jet quantities: 

[ To exploit the increased experimental precision reliable theoretical predictions 

in next-to-leading order (NLO) of QCD perturbation theory are required. In 
this proceedings contribution I give an overview of recent developments in NLO 
calculations for deeply-inelastic processes. The main improvement during the 
last two to three years was that universal Monte Carlo programs have become 
available which permit the numerical calculation of any (2-f l)-jet-like infrared- 
safe observable in NLO. 

• Matching of DIS and photoproduction: 

For photoproduction (Q^ « 0) and deeply inelastic scattering (DIS, ^ Aqcd) 
it is well known how to calculate cross sections systematically in perturbation 
theory. Recently, a formalism ^ has been developed which permits calculations 
in the transition region ^ Aqcd- 
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• One-particle-inclusive processes: 

The measurements of momentum fraction distribution mentioned above show a 
very good agreement of experimental data and theoretical predictions in next-to- 
leading order for moderately large and Xp. However, the theoretical prediction 
breaks down both for small and small Xp. 

For lack of space, I had to leave out many, if not most, interesting topics. For 
more detailed information, I would like to refer the interested reader to the proceedings 
of the DIS 98 workshop in Brussels 0. 

2. NLO calculations for jet quantities 

One of the basic problems of perturbative QCD calculations is that experimentally 
hadrons are observed in the final state while theoretical calculations yield results for 
partons. Moreover, not all observables can be calculated in perturbation theory in a 
meaningful way. In principle, there are two possibilities: 

• Infrared- safe observables, which are constructed such that all soft and coUinear 
singularities cancel among real and virtual corrections or can be absorbed into 
redefined parton densities. The same observable is then evaluated both for 
parton final states (theory prediction) and hadron final states (experimental 
data), possibly after the experimental data have been corrected for systematic 
errors. 

• Alternatively, additional non-perturbative objects can be introduced, for instance 
fragmentation functions, which allow for a study of one-particle-inclusive 
processes. Fragmentation functions have to be measured and parametrized 
experimentally, and may serve to hide final-state coUinear singularities which 
do not cancel because of integrations over restricted phase space regions. 

In this section I consider the first possibility; one-particle-inclusive processes will be 
treated in Section ^. 

In QCD perturbation theory, expectation values for parton observables are 
calculated as a phase space integral of a product of a differential parton cross section 
c'-"'' (pii ■ ■ ■ ,Pn) for n-parton final states and an observable O^"^ {pi, . . . ,Pn)'- 

{0)=Y, f dPS(") (pi , . . . , p„) (Pi , . . . , P„) . (1) 
In next-to-leading-order calculations, there are three contributions to be included: 

{O) - aBorn O^-^-^^ + ^virtual O^""'' + fT^eal O^"' . (2) 

Here (TBorn is the lowest order cross section, (Tvirtuai are the virtual and fircai are the 
real corrections. If the Born term has n — 1 final-state partons, (ivirtuai will also have 
n — 1 and (Ticai will have n final-state partons. Infrared singularities arise in (Tvirtuai in 
the loop integrations and in (Tioai in the phase space integration over dPS''"^ 

As already mentioned in the introduction, theoretical predictions for partons in 
the final state are infrared-finite only for a special class of observables. The technical 
requirement for infrared-safe observables is that they behave well under soft and 
coUinear limits: 

(pi, . . . . . . ,Pn) O^"-^) (Pl, . . . . . . ,p„) , 

(pi, ...,p„...,pj,...,pn) > 0^""^) (pi, ...,p„...,pj,... ,Pn,Pi + 

PiWpj 
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Momenta denoted by p are to be omitted. 

The main technical problem is the extraction of the infrared singularities from 
the real corrections. It turns out that this can be done in an observable-independent 
way, such that it is possible to build Monte-Carlo programs which are able to integrate 
arbitrary infrared-safe observables. This can be done because the structure of QCD 
cross sections in kinematical limits is known: the factorization theorems of QCD 
state that the structure of the parton cross section CTrcai for coUinear and soft limits is 
of the form of a product of a singular kernel K and the Born cross section ct'"^^-': 

fjf") _^ Ka^'"^\ (4) 

soft/collincar 

The product of CTicai and O*^"^ thus behaves in a simple way: the cross section goes 
over into a kernel K and the Born cross section, and the observable approaches the 
corresponding observable for Born term kinematics. The kernel K is independent of 
the phase space variables of the {n — l)-particle phase space, and thus the phase space 
integration over the corresponding variables can be performed analytically. 



2.1. Calculations 

Particularly interesting for phenomenological applications are processes with 2-f 1 jets 
in the final state (which means that 2 jets are produced from the hard scattering cross 
section, plus the remnant jet of the incident photon). First of all, in leading order of 
QCD perturbation theory these processes are of 0{as), and are thus suitable for a 
measurement of the strong coupling constant as. Moreover, the gluon density enters 
in leading order in the so-called boson-gluon-fusion process. Therefore, this process 
can also be used to measure the gluon density |l^ . 

By now there are several calculations for (2-t-l)-jet processes available with 
corresponding weighted Monte-Carlo programs: 

• PRO JET jl4j: The jet definition is restricted to the modified JADE jet clustering 
scheme; the program is based on the calculation published in Rcfs. |l^, p^ . 

• DISJET |l^: Again the jet definition is restricted to the modified JADE scheme; 
the program is based on the calculation in Rcfs. |2C| |. 

• MEPJET [|T]: This is a program for the calculation of arbitrary observables 
which uses the phase-space-slicing method. The corresponding calculation 
uses the Giele-Glover formalism ||23j for the analytical calculation of the IR- 
singular integrals of the real corrections, and the crossing- function technique 

to handle initial-state singularities. The latter requires the calculation of "crossing 
functions" for each set of parton densities. 

• DISENT |2^: This program is based on the subtraction method. The subtraction 
term is defined by means of the dipole formalism| ^ . 

• DISASTER++ nil: This is a C++ class libraryfj] . The subtraction method is 
employed, and the construction of the subtraction term resembles the method 

§ The subtraction term is written as a sum over dipoles (an "emitter" formed from two of the 
original partons and a "spectator" parton). Besides the factorization theorems of perturbative QCD, 
the main ingredient is an exact factorization formula for the three-particle phase space, which allows 
for a smooth mapping of an arbitrary 3-parton configuration onto the various singular contributions. 
II The acronym stands for "Deeply Inelastic Scattering: All Subtractions Through Evaluated 
Residues". Most of the program is written in C++. A FORTRAN interface is available; thus there 
is no problem to interface the class library to existing FORTRAN code. 
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of Ref. [£9j , i.e. it is obtained by the evaluation of the residues of the cross section 
in the soft and coUinear hmits. Double counting of soft and collinear singularities 
is avoided by means of a general partial fractions method. 

• JetViP This program implements the calculation of |3l], which extends the 
previous calculations into the photoproduction limit Q'^ ^ 0. The calculation 
has been done by means of the phase space slicing method. Up to now, the 
polarization of the virtual photon is restricted to be longitudinal or transverse. 

The two basic approaches which are employed to extract the infrared singularities from 
the real corrections are the phase-space-slicing method and the subtraction method. 

• The phase-space-slicing method splits up the full parton phase space into two 
regions: a region R where all partons can be resolved, and a region U where 
two or more partons are unresolved. This split is usually achieved by means 
of a technical cut parameter Smm- Two partons with momenta pi and p2 are 
unresolved if their invariant mass 2pip2 is smaller than Smin and resolved if it 
is larger. The integration over the resolved region R can be performed safely 
by Monte Carlo integration, because all infrared singularities are cut out by the 
phase space cut. The integration over the unresolved region U is divergent and 
cannot be performed numerically, but because of the constraint 2pip2 < Smin the 
cross section factorizes (see Eq. ^ in the limit Snun 0. This contribution is 
approximated by this limit. The integration over the singular region can be done 
analytically, and the divergent parts can be extracted. In the limit of Smin 
the sum of the two integrals over R and U should approach the integral over 
the full phase space. It should be kept in mind that this convergence has to be 
checked explicitly by varying s,„in and looking for a plateau in this variable. 

• A calculation using the subtraction method defines a subtraction term S which 
makes the integral J dPS (cr^"^ O*^"^ — S) finite. The original integral is, as an 
exact identity, rewritten as 

J dPSd*") = J dPS ((t(") - S*) dPSS*. (5) 

The first integral can be done by a Monte Carlo integration. For the term S, the 
factorization from Eq. ^ holds exactly. As for the phase-space-slicing method, the 
second term is integrated analytically. No technical cut-off has to be introduced^. 

Both methods have their merits and their drawbacks. The phase-space-slicing method 
is technically simple and can be easily implemented once the matrix elements for the 
real and virtual corrections are known. The main problem is the residual dependence 
on the technical cut Smin- The independence of numerical results from variations of this 
cut has to be checked; moreover, the integration over the region R mentioned above 
requires very high statistics, because the integration region is close to the singular limit. 
The subtraction method does not require a technical cut, but the construction of the 
subtraction term S is usually quite involved. If this can be afforded, the subtraction 
method is the method of choice. 

^ This is, strictly speaking, not correct. A dimensionless cut t(.ut of the order of 10"^" to 10~^^ 
is used to avoid phase space regions where the subtraction no longer works because of the finite 
precision of floating point numbers. 



Jets and fragmentation 



5 



2.2. Program comparisons 

It is interesting to compare the available universal Monte Carlo programs numerically 
to check whether all available calculations are consistent. Experimental papers usually 
contain statements that the programs "agree on the one per cent level". A closer 
investigation, however, reveals that a statement of this kind is not correct. The three 
programs MEPJET, DISENT and DISASTER++ have been compared in Ref. ^ for the 
modified JADE jet clustering algorithm in the E-scheme for several choices of physical 
and unphysical parton densities^ The result is that DISENT . 1 and DISASTER++ 1 . 
agree well, with discrepancies of the MEPJET results. Presently this is studied in the 
framework of the HERA Monte Carlo workshop]^ at DESY. There does not yet exist a 
systematic comparison of the JetVip program with MEPJET, DISENT and DISASTER++. 



3. Matching of DIS and photoproduction 

For large photon virtuality Q^, the coupling of the exchanged virtual photon in a 
lepton-nuclcon scattering process is exclusively pointlike. Extending this kind of 
calculation down to « leads to the problem that the photon propagator diverges. 
Instead, it is possible to calculate the scattering process for the scattering of a quasi- 
real photon and a nucleon, where the flux of quasi-real photons is described by a 
Weizsacker- Williams approximation. For small Q^, in addition to the cross section 
contribution from the pointlike coupling, a resolved contribution has to be added, 
because the quasi-real photon may fluctuate into a hadronic state, which in turn 
interacts strongly with the incident nucleon. This process is modelled by means of 
parton densities fi/^ of the virtual photon. The assumption of a photon structure is 
also required in order to treat coUinear singularities arising from the splitting of the real 
photon via its pointlike coupling into a coUinear quark-antiquark pair. This coUinear 
singularity does not cancel against the virtual corrections, but is absorbed into the 
fi/ry. Typically, the infrared singularities are regularized by dimensional regularization; 
the singularities then show up as poles in e, where the space-time dimension is set to 
d = 4-2e. 

This type of calculation has recently been extended to the case of exchanged 
photons with moderate in Refs. |^. Here, because is finite, strictly speaking 
there is no coUinear singularity, and therefore no poles in e related to the photon 
splitting arise. However, the integral over the phase space of the quark-antiquark pair 
yields a logarithm in Q^. Because this logarithm may be large, and can therefore spoil 
perturbation theory, it has to be resummed. This is done by absorbing it into the 
redefined parton densities of the photon. The corresponding renormalization group 
equation then takes care of the resummation. 

Depending on the factorization scale for the virtual photon, the resolved 
contribution can be surprisingly large even for fairly large Q^, compared with 
the "standard" DIS calculation for the pointlike coupling. This seems to be in 
contradiction with the statement that the resolved contribution should die out for 
increasing Q^. There are two reasons for this: (a) The choice of the factorization 
scale fj, for the resolved photon dictates the size of the resolved contribution. The 

+ By "unphysical" I mean parton densities of the form q(x) = (1 — x)" and g{x) = (1 — x)" , where 
a is some power. These are introduced to have a more stringent test on the hard scattering matrix 
plpments 
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parton density of the virtual photon is exactly zero for fi — Q. The factorization 
scale employed in Refs. ^ ^ is \/Q^ + (Et is the transverse energy of the 
produced jets), which makes sure that even at large there is always a resolved 
contribution, (b) In the full NLO calculation, there are four different matrix elements 
which contribute: the direct process in LO and NLO (with the dangerous logarithm 
in subtracted), and the resolved process in LO and NLO. It is expected that the 
sum of the first three processes reproduce the result for the standard calculations, 
and this is indeed the case: the logarithm that has been subtracted for the direct 
coupling is added up again via the parton density of the photon in the LO resolved 
contribution. The difference comes from the resolved contribution in NLO: the 
corresponding parton subprocess, which is of O (a^), is one order in Us higher than 
the "standard" calculation. Thus, this contribution could be considered as part of the 
NNLO correction to the Born term for (2+l)-jet production. Differences between the 
two approaches are therefore expected. 

4. One-particle-inclusive processes 

The comparison of Xp-distribution^ for charged particle production from experimental 
data 1^, and the NLO program CYCLOPS ||3|, Q leads to severe discrepancies for 
small values of or small Xp. For large and large Xp, data and theory agree 
nicely. Where does this discrepancy come from? 

The theoretical prediction is made in the fragmentation function picture: the cross 
section for the inclusive production of charged particles is obtained by a convolution 
of the hard scattering cross section calculated in perturbative QCD and fragmentation 
functions which have been obtained from fits to e"'"e~ data. Fragmentation functions 
depend on the momentum fraction z of the parent parton carried by the observed 
particle. An assumption in this picture is that the mass of the observed particle 
can be neglected relative to any other scale of the process, in particular relative to its 
momentum. The variable z can thus be defined either by means of fractions of energies 
or fractions of momenta. In the real world, the observed particle has a mass, and this 
gives, thus, rise to an uncertainty in the theoretical description. It is clear that mass 
effects will be important if Xp — 0(2771^/(5), being a typical hadronic mass. It 
turns out that excluding data points with a value of Xp close to or smaller than this 
leads to a good agreement between data and theory. A different argument in terms 
of rapidities of partons and observed particles has been given in Ref. jsj] . During the 
Durham workshop, a quantitative estimate of power corrections ~ to the fixed- 

order NLO prediction has been made. Y. Dokshitser and B. Webber proposed a factor 
l/(l + 4/i^/ {xpQ)'^), depending on a mass parameter /x, to be multiphed with the NLO 
cross section; this factor together with a fit of /x is able to describe the experimental 
data fairly well (see the contribution to these proceedings by P. Dixon, D. Kant, and 
G. Thompson). 

5. Summary 

I have discussed three topics related to hadronic final states at HERA. For the 
basic processes, theoretical predictions are available in next-to-leading-order accuracy. 

tt The variable Xp is defined to be the fraction 2E/Q, where E is the energy of an observed particle 
in the current hemisphere of the Breit frame. 
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Independent calculations permit the comparison of results, and a few problems with 
Monte Carlo programs have already been fixed. What is still missing are calculations 
for W and Z-exchange in the subtraction formalism for jet cross sections. This is 
likely to become available in the near future. Moreover, a calculation for transverse 
momentum spectra of charged particles has not yet been done. The calculation for the 
transition region of DIS and photoproduction fills a gap in the theoretical description 
of lepton-nucleon scattering. However, it is not yet clear whether the parton densities 
for virtual photons are process- independent beyond NLO, such that they can be 
measured in one process and used for predictions in a different one. The necessity 
to introduce a power correction term for one-particle-inclusive distributions already 
at fairly large values of shows that the calculation of fixed-order QCD corrections 
is not sufficient for a good description of experimental data. Unfortunately, a power 
correction term introduces an additional mass parameter, which cannot be calculated 
from first principles. 
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